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Abstract. We describe the analysis of longitudinal air shower profiles as measured by the fluorescence detectors of the Pierre 
Auger Observatory and present the measurement of the depth of maximum of extensive air showers, X^ax, with energies > 
10 eV. The measured energy evolution of the average of X^ax and its fluctuations, RMS(X[nax)> are compared to air shower 
simulations for different primary particles. 
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INTRODUCTION 

The determination of the chemical composition of ultra- 
high energy cosmic rays is essential to understand the 
origin of cosmic rays and to interpret the features ob- 
served in the ultra-high energy cosmic ray flux. For in- 
stance, the observed hardening of the cosmic ray en- 
ergy spectrum at energies between 10^^ eV and lO'^ eV, 
known as the 'ankle', might either be a signature of the 
transition from galactic to extragalactic cosmic rays or a 
distortion of a proton-dominated extragalactic spectrum 
due to energy losses [1]. Moreover, the flux suppres- 
sion observed above 4-lO'^eV [2] could be either due 
to propagation effects [3] (photopion production of pri- 
mary protons or photonuclear reactions of primary nu- 
clei) or a signature of the maximum injection energy of 
the sources [4]. 

There are several experimental methods to estimate 
the primary composition from cosmic ray induced air 
showers. Within the Pierre Auger Observatory (see [5] 
and [6]) the observation of the longitudinal shower de- 
velopment with fluorescence detectors allows to measure 
the depth of the maximum of the shower evolution, X^ax, 
which is sensitive to the primary mass^ 

With the generalization of Heitler's model of electron- 
photon cascades to hadron-induced showers [8] and the 
superposition assumption for nuclear primaries of mass 
A, the average depth of the shower maximum, {Xmax), at 
a given energy E is expected to follow 



where (InA) is the average of the logarithm of the pri- 
mary masses. The coefficients a and j3 depend on the 
nature of hadronic interactions, most notably on the mul- 
tiplicity, elasticity and cross-section in ultra-high en- 
ergy collisions of hadrons with air, see e.g. [9]. The 
change of (Xmax) pei" decade of energy is called elon- 
gation rate [10], Diq — and it is sensitive to 
changes in composition with energy. A complementary 
composition-dependent observable is the magnitude of 
the shower-to-shower fluctuations of the depth of max- 
imum, RMS(X,Tiax), which is expected to decrease with 
the number of primary nucleons A (though not as fast 
as 1/VA [11]) and to increase with the interaction 
length of the primary particle. In case of a mixed com- 
position, the full width of the X^^x distribution follows 
from the shower-to-shower fluctuations of the individ- 
ual mass groups and their separation in (Xmax) [12]. For 
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' For other methods based on the surface detector see e.g. [7]. 



FIGURE 1. RMS(Xmax) from different hadronic interaction 
models [23] and a two-component p/Fe composition model 
(£=10"*eV). 




FIGURE 2. Left: Illustration of the field of view bias. Right: Dependence of (Xmax) on the field of view boundary close to ground 
for data and MC {10^^ '^ eV < £ <10^^ ^ eV). 



a simple two-component composition of primaries with 
masses A 1 and A2 and abundances /i and /i = (1 — /i), 
RMS(Xmax) is given by 

RMS(X„ax) = (/iVi +/2V2 + /i/2Ai)^ (2) 

where Vi = RMS(Xn,ax)^ and Ax = (^max)i - (^max)2- 
As can be seen in Fig. 1, a proton/iron mixture gives rise 
to a broad maximum in RMS(X,^ax) for 0.3 < /p < 1 and 
a rapid decrease of the width towards /p = 0. 

DATA ANALYSIS 

The Xjnax results which were discussed at this workshop 
are based on [13] and use air shower data recorded be- 
tween December 2004 and March 2009. Only events 
detected in hybrid mode [14] are considered, i.e. the 
shower development must have been measured by the 
fluorescence detector (FD), and at least one coincident 
surface detector station is required to provide a ground- 
level time. Using the time constraint from the surface 
detector, the shower geometry can be determined with 
an angular uncertainty of 0.6° [15]. The longitudinal 
profile of the energy deposit is reconstructed [16] from 
the light recorded by the FD using the fluorescence and 
Cherenkov yields and lateral distributions from [17]. 
With the help of data from atmospheric monitoring de- 
vices [18] the light collected by the telescopes is cor- 
rected for the attenuation between the shower and the 
detector and the longitudinal shower profile is recon- 
structed as a function of atmospheric depth. Zmax is de- 
termined by fitting the reconstructed longitudinal profile 
with a Gaisser-Hillas function [19]. 



To assure a good Zmax resolution, the following qual- 
ity cuts are applied: The impact of varying atmospheric 
conditions on the Xmax measurement is minimized by re- 
jecting time periods with cloud coverage and by requir- 
ing reliable measurements of the vertical optical depth 
of aerosols. Profiles that are distorted by residual cloud 
contamination are rejected by a loose cut on the qual- 
ity of the profile fit ix^fNdf<2.5). We take into account 
events only with energies above lO'^ eV where the prob- 
ability for at least one triggered surface detector station 
is 100%, irrespective of the mass of the primary parti- 
cle [20] . The geometrical reconstruction of showers with 
a large apparent angular speed of the image in the tele- 
scope is susceptible to uncertainties in the time synchro- 
nization between the fluorescence and surface detector. 
Therefore, events with a light emission angle towards the 
FD that is smaUer than 20° are rejected. This cut also re- 
moves events with a large fraction of Cherenkov Ught. 
The energy and shower maximum can be reliably mea- 
sured only if Xmax is in the field of view (FOV) of the 
telescopes (covering 1.5° to 30° in elevation). Events 
for which only the rising or falUng edge of the pro- 
file is detected are not used. Moreover, we calculate the 
expected statistical uncertainty of the reconstruction of 
^max for each event, based on the shower geometry and 
atmospheric conditions, and require it to be better than 
40g/cm2. 

The latter two selection criteria may cause a system- 
atic under- sampling of the tails of the Xmax distribution, 
since showers developing very deep or shallow in the at- 
mosphere might be rejected from the data sample (see 
illustration in the left panel of Fig. 2). To avoid a corre- 
sponding bias in the measured {Xmax) and RMS(Xmax) 
we apply fiducial volume cuts on the viewable Xmax 




FIGURE 3. (Xtnax) and RMS(Xn,ax) compared with air shower simulations [22] using different hadronic interaction models [23]. 



range. For this purpose the effective upper and lower field 
of view Umits are calculated for each event and (^max) is 
measured as a function of these Umits. An example of 
the (Xmax) dependence on the lower FOV limit is shown 
in Fig. 2 for data and simulated events. As can be seen, 
the (Xniax) is asymptotically unbiased for events with a 
sufficiently deep FOV limit, but it is systematically too 
shallow when the FOV boundary starts cutting into the 
tails of the X^ax distribution. Obviously, the unbiased re- 
gion depends on the distribution itself and can thus not 
be determined by simulations. Instead, we fit the data 
with the mean of a one-sided truncated normal distribu- 
tion (shown as solid lines in Fig. 2) and reject all events 
that have a FOV limit for which the measured (Xmax) de- 
parts by more than 5 g/cm^ from its asymptotic value. 

After all cuts, 3754 events are selected for the X^m. 
analysis. The X^ax resolution as a function of energy for 
these events is estimated using a detailed simulation of 
the FD and the atmosphere. The resolution, defined by 
the full standard deviation, is at the 20 g/cm^ level above 
a few EeV. The difference between the reconstructed 
^max values in events that had a sufficiently high energy 
to be detected independently by two or more FD stations 
is used to cross-check these findings and as it was shown 
in [13], the simulations reproduce the data well. 

RESULTS 

The measured (Xmax) and R]V[S(Xmax) are measured in 
energy bins of AlgE = 0.1 below 10 EeV and AlgE = 
0.2 above that energy. The last bin starts at lO'^'^eV, 
integrating up to the highest energy event (E = (59 ± 
8) EeV). The systematic uncertainty of the FD energy 
scale is 22% [24]. Uncertainties of the calibration, at- 
mospheric conditions, reconstruction and event selection 
give rise to a systematic uncertainty of <13g/cm^ for 
(^max) and <6g/cm^ for the RMS. The results were 
found to be independent of zenith angle, time periods and 



FD stations within the quoted uncertainties. 

The measured values are displayed in Fig. 3. A fit of 
(-^max) data with a constant elongation rate does not de- 
scribe our data (;if^/Ndf=34.9/ll), but using two slopes 
yields a satisfactory fit (;t^/Ndf=9.7/9) with an elonga- 
tion rate of (106^^^) g/cm^/decade below loi8-24±o.05 
and (24±3) g/cm^/decade above this energy. If the prop- 
erties of hadronic interactions do not change significantly 
over less than two orders of magnitude in primary energy 
(< factor 10 in center of mass energy), this change of 
ADio =(82^2i) g/cm^/decade would imply a change in 
the energy dependence of the composition around the an- 
kle, supporting the hypothesis of a transition from giilac- 
tic to extragalactic cosmic rays in this region. 

The shower-to-shower fluctuations, RMS(Xniax), are 
obtained by subtracting the detector resolution in quadra- 
ture from the width of the observed X^ax distributions 
resulting in a correction of <6g/cm^. As can be seen in 
the right panel of Fig. 3, we observe a decrease in the 
fluctuations with energy from about 55 to 26 g/cm^ as 
the energy increases. Assuming again that the hadronic 
interaction properties do not change much within the ob- 
served energy range, these decreasing fluctuations are an 
independent signature of an increasing average mass of 
the primary particles. 

For the interpretation of the absolute values of (Xmax) 
and RMS(Xmax) a comparison to air shower simulations 
is needed. As can be seen in Fig. 3, there are considerable 
differences between the results of calculations using dif- 
ferent hadronic interaction models. These differences are 
not necessarily exhaustive, since the hadronic interaction 
models do not cover the full range of possible extrapo- 
lations of low energy accelerator data. If, however, taken 
at face value, the comparison of the data and simulations 
leads to the same conclusions as above, namely a gradual 
increase of the average mass of cosmic rays with energy 
up to 59 EeV. 

It is illustrative to compare the data with predictions 
for a simple two-component proton/iron model using 




FIGURE 4. (Xmax) and RMS(Xmax) compared to QgsjetII predictions for a two-component iron/proton composition. Data 
points show extreme values obtained by shifting the measurements by ± one sigma of the systematic uncertainties. 



the QgsjetII hadronic interaction model. (Xmax) and 
RMS(Xmax) simulations for different proton fractions / 
are shown in Fig. 4. As can be seen, the (Xmax) values 
change linearly with / as expected from Eq. (1), whereas 
the width of the Xmax distribution is very similar for pro- 
ton fractions f ^ 03. To visualize the systematic un- 
certainties of the data, in this figure we shifted the de- 
fault results from Fig. 3 by their systematic uncertainties. 
Note that the systematics on RMS(Xmax) are dominated 
by the event selection, whereas the systematics on (Xmax) 
are mainly due to reconstruction uncertainties and atmo- 
spheric effects. Therefore sys((Xmax)) and sys(RMS) are 
uncorrelated and can be shifted independently. Within 
this simplistic two-component model, the data is com- 
patible with a light or mixed composition at low energies. 
At high energies, a heavy composition would result, but 
the (Xmax) would indicate a larger proton fraction than 
RMS(Xniax)- At high energies, this model corresponds to 
a heavy composition, however, the (Xmax) would indicate 
a larger proton fraction than RMS (Xmax) • 
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